Magnetic Gaps related to Spin Glass Order in Fermionic Systems 
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We provide evidence for spin glass related magnetic gaps in the fermionic density of states 
below the freezing temperature. Model calculations are presented and proposed to be relevant 
for explaining resistivity measurements which observe a crossover from variable-range- to activated 
behavior. The magnetic field dependence of a hardgap and the low temperature decay of the density 
of states are given. In models with fermion transport a new metal-insulator transition is predicted 
to occur due to the spin-glass gap, anteceding the spin glass to quantum paramagnet transition at 
smaller spin density. Important fluctuation effects due to flnite range frustrated interactions are 
estimated and discussed. 
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Magnetic hard gaps were observed in a large number 
of systems such as In-doped CdMnTe , amorphous 
GeCr-films boron- and arsenic-doped siUcon [^-0, 
in amorphous S'li-^jMn^j-films and irradiated poly- 
mers [p|JlO|] by manifestation of a thermal crossover from 
variable range hopping- or Mott-behavior to activated 
ea::p(-Bi//T)-behavior in the low T resistivity. Reduction 
of the activation energy associated with the gap under 
an increased external magnetic field was held responsi- 
ble for both the generally observed feature of a negative 
magnetoresistance and the magnetic origin of the gap. 
In all experiments the crossover in resistivity behavior 
was seen at temperatures close to the corresponding spin 
glass freezing temperature. An existing theory involving 

which explained the 



localized magnetic polarons 1 1 1 



gradual hardening of Coulomb gaps was consid- 

ered to be a possible and appropriate explanation of some 
of the experiments [^] . For other experiments P,p|,p|pl)| 
spin glass order was suggested to be at the origin of the 
observed activated behavior. The random position of the 
implanted ions in these materials leads to a broad distri- 
bution of exchange integrals which favors random freez- 
ing of the magnetic moments. 

Fermionic spin glasses are marked by the fascinating 
combination of a complex energy landscape with phe- 
nomena well known from strongly correlated systems 
like quantum phase transitions, non Fermi liquid behav- 
ior, phase separation, and interaction induced metal- 
insulator transitions. The paramagnet to spin glass quan- 
tum phase transition occuring in these systems was stud- 
ied in [T^-[l7|] , while strong correlation effects deep inside 
the spin glass phase were discarded. Recently the proper- 
ties of highly correlated metals and interaction-induced 
insulators have attracted much interest |Q . 
The goal of this Letter is to report theoretical results, ob- 
tained by replica-methods and TAP-equations, for com- 
parable phenomena in several different fermionic spin 
glasses, which develop either a hard gap or a pseudogap 
in the fermionic density of states at zero temperature. 



depending in size and shape on the spin glass order pa- 
rameter(s). In any case the manifestation of spin glass 
order in this fermionic property represents a surprisingly 
large effect. 

The magnetic field induced reduction of the gap, im- 
plying a negative magnetoresistance, and its stability 
against fermion hopping transport are demonstrated. 
Since the system stays insulating due to spin glass or- 
der it may be characterized as a SG (spin glass) based 
insulator. The existence of a SG based metal insulator 
transition is conjectured from the combination of our new 
calculations with results obtained in the vicinity of the 
quantum paramagnet to spin glass transition in metallic 
systems . The low energy behavior of the Green func- 
tion in the SG based insulating phase is derived and an 
estimation for the location of the abovementioned tran- 
sition is given. We discuss implications of our results for 
the interpretation of the experiments [p]-p^ . 
The fermionic Ising spin glass ISGf is defined by the 
Hamiltonian Ti, = Jij<^i<^j — fJ-^ni with spins 

a ~ nij — n^, the fermion number operator n — ni; + n^^, 
chemical potential fj, , and Gaussian-distributed Jij with 
variance J^. It is the simplest model that takes charge- 
spin couplings and corresponding fluctuations into ac- 
count and, by virtue of its exact solvability in infinite 
space dimensions and for T — > provides solid theoreti- 
cal background to be compared with related experimental 
situations. Due to the possibility of double and zero oc- 
cupancy of sites the freezing temperature is slightly lower 
than the one of its standard counterpart (SK-model) 
with one spin at each site. Thermally activated hop- 
ping of the localized fermions is effectively allowed for 
by working in the grand canonical ensemble with statis- 
tically fluctuating particle numbers at each site. While 
spin correlation functions of the model are static, corre- 
lators with an odd number of equal time fermion oper- 
ators display the rich quantum-dynamical structure of 
an interacting fermion system. The fermionic path in- 
tegral technique provides the tool for the calculation of 
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quantities such as density of states, (dc or ac) conduc- 
tivities et cetera for all kinds of spin glass models with 
charge degrees of freedom. We consider important the 
question, whether such a purely magnetic model, whose 
Hamiltonian consists just of a frustrated spin interaction, 
can have specific and strong effects on charge correlations 
and fluctuations, and whether a continuous distribution 
of exchange energies can cause a rigid gap of finite width. 
The density of states /Co-(e) — ~^IniG^{e), calculated by 
means of the replica-method from the fermion Green's 
function of the ISGf , reads 



ch{P^j.) +chij3{e + fi)) 



dz 



-^-(3{e + fi + aH{z)r)/{2J^x) 



(1) 



where H{z) = H+Jy^z,H, q, x denote effective field, ex- 
ternal magnetic field, spin glass order parameter, and lin- 
ear susceptibility. Effects from replica symmetry break- 
ing (RSB), involving oo-many order parameters qi and 
variables Zi, and from finite-range interactions Jij are 
not yet included here, but are considered below. The 
existence of a hardgap of size 2Eg{H) is nevertheless an 
unexpected and surprisingly large effect of spin glass or- 
der; it is best appreciated in our result for the density of 
states result at T = 0, which is given by 

pAE)\^^, = e , (2) 



27tJ 

with E := e + /i and the field-dependent gap energy 
Eg{H) = J^y2/^exp{-H^/{2J^)). Eq.(|) is valid within 
the regime \^\ < Eg/2, which corresponds to a half filled 
system at T = with filling factor ly = 1 + {1 — q — 
Tx)tanh{Pfi). In terms of Eg{H) the following low tem- 
perature expansions (for linear susceptibility and order 
parameter), required for the exact evaluation of the DoS 
formula (|l|) at low T, were obtained as q = 1 — Eg{H)T + 
0{T^),X^ Eg[H) + \El{H){l-H^)T+0{T^). The last 
equation allows one to express the zero field gapwidth as 
Eg = a{^)xT^ completely in terms of experimentally ac- 
cessible quantities, where a increases with the chemical 
potential from 2.1833 at ^ = to 2.4075 at = J/%/27r. 
These results show that i) the density of states is zero at 
T = in the finite interval given by |e -|- /Lt| < Eg{H) and 
ii) the gapwidth in a magnetic field is independent of spin 
orientation (in contrast to the local-limit Hubbard gap 
in an external field) and shrinks as the field is increased. 
This will be the source of the negative magnetoresistance 
in the extended Ising spin glass model with charge trans- 
port as discussed below. The physical consequences of 
these spin glass related properties are hence in agreement 
with the experimental observations of crossover behavior 
in the low T resistivities mentioned above. For single 



fermion energies E = e + smaller than the gap energy 
Eg{H) and \fi\ < \Eg{H) the density of states decays to 
zero exponentially as given by 
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2 J cos{^- 



{l3Eg{H))-^ [chiPtx) + ch{l3E)] 



g Eg{H) 



(3) 




FIG. 1. Single particle density of states (DoS) for the 
fermionic Ising spin glass (ISGf) versus energy and magnetic 
field in units of J. 

The infinite range insulating spin glass model allows 
to evaluate effects of replica symmetry breaking (RSB). 
Let us discuss the {T = 0, H = 0) result for the fermion 
Green's function 

G^'iE) = E / ^(^) • (4) 

A=±l"'" E+iO-X[h+x] 

The weight P{H) depends on the set of Parisi order pa- 
rameters and simplifies to e~~ under replica sym- 
metry reproducing thus Eq.^ together with x = X- 
Hence the one particle excitations reveal a gap if the 
nonequilibrium susceptibility Xj which enters as a self 
energy and represents the effect of the Onsager reaction 
field, is finite. The influence of the randomly distributed 
local fields on the other hand is contained in the weight 
P{H). For the infinite-range model, the RSB-solution 
shows that x vanishes like T for T — s- 0. Our explicit cal- 
culations show that in this case the hardgap turns into a 
pseudogap (details will be published elsewhere). Finite 
range interactions however lead to corrections that ren- 
der X finite [|l|, which indicates that fluctuations will 
tend to reinstall the gap. Whereas a controlled theory 
for self energy corrections to Eq.(^ does not yet exist, 
an estimation of these small fluctuation contributions to 
the density of states is given below. 
The gap also turns out to be robust against an addi- 
tional fermion hopping term until the bandwidth exceeds 
a critical value. The gap and related properties of this ex- 
tended itinerant spin glass model, whose magnetic phase 
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diagram and -transitions have previously been analyzed 
1^,0, are also discussed below. 

Coulomb interaction effects are of diverse nature: the 
long range part will still tend to depress the (remain- 
ing) DoS near the Fermi level and thus stay particularly 
relevant when either spin glass order is weak or absent, 
or (in the case of fully developed spin glass order) when 
the Fermi energy lies close to the gap edges. The Hub- 
bard coupling U leads to a shift of the chemical potential 
{n fj,— ^v), and also to a shift H H' = H+^m{H) 
of the applied magnetic field. 

Corrections associated with a finite interaction range are 
estimated by considering TAP-equations for the ISGf. 
Similarly to the spin case an Onsager reaction term 
— •A^o'iXjj™! (coordination number z) must be 

subtracted from the local field on the level of frozen in 
magnetizations, but to account for spin fluctuations a self 
energy ^ Jfj'^iXjj'^i must be added again. Thus the 

polarization (JiJijXjj of surrounding spins lowers the en- 
ergy of magnetic sites by Eg = Ej=i Xiil-^ijO^la-u- As 
this gap energy originates in the formation of polarized 
spin clusters it differs from the Hubbard gap width U 
caused by onsite repulsion. Now we argue that for a fi- 
nite interaction range the local susceptibility behaves like 
X = const. + 0{T) instead of the x ~ found by solving 
the TAP-equations for the infinite-range SK-model . 
The latter is valid for a system stuck in a fixed free en- 
ergy valley with infinitely high barriers and is intimately 
related to a pscudogap P{hioc) ^ h in the distribution of 
local fields, which is absent in finite dimensional systems 
pl| . Thus fiuctuation effects in systems with finite-range 
frustrated interactions must be expected (despite the not 
yet available finite-range spin glass theory) to favor the 
hardgap as predicted by the replica symmetric solution 
eq.(y) of the infinite range model. Hence we conjecture 
p{e) ~ [9(e - xZ]j=i(-^y )^)]au ^ and therefore along 
the lines of Q anc exp(— i?/T~) for bonds with a 
Gaussian distribution. A sufficiently large (magnetic) co- 
ordination number z or a long-ranged RKKY-interaction 
leads to an activated behavior of the conductivity, other- 
wise nonuniversal temperature exponents may occur. A 
bimodal ± J - distribution on the other hand is expected 
to result in a true hardgap. 

Multi-valley correlations of statistical DoS - fluctuations 

.(fc) - 



phase at II = 0.9, but when analyzing the stability of the 
magnetic solution according to the scheme of de Almeida 
and Thouless |^ we found two negative eigenvalues of 
the Hessian matrix instead of the expected one negative 
value indicating just the instability towards RSB. As de- 
scribed in 1^ for spin glass problems, one has to choose 
the stable solution with the lowest free energy, hence the 
system undergoes a first order transition from half filling 
to the completely filled paramagnetic state at = ■^Eg. 
The electronic kinetic energy is properly embedded in an 
itinerant /S'G'/-model by adding the random hopping 

Hamiltonian Ht = J2ij-a ^ij'^l a^j,'^- The real-symmetric 
tij are Gaussian-distributed and chosen to describe a 
semi-circular band of width Eq. The 4-fermion term 
which results from -averages is decoupled by means of 
a quaternionic matrix field i? in a way familiar from An- 
derson localization theory. The saddle point value ra{ei) 
of this field obeys the selfconsistency equation r^iei) = 

'-f cAei) = -'-BEalodre^;'^ < T.c,^.(r)ct ,,(0) >, 
which becomes exact for infinite space dimensions. The 
expectation value employs the quantum-dynamical ac- 

tion A = Jo^ Jo^dTdT'{j:^^^cl,{TmT ~ T'){dr + /i) + 
ira[T,T')]CaAr) + Ea,b'5"(T)5''(T')Q'"^(T',T)) COutaiu- 

ing the disorder averaged spin interaction. In contrast 
to the quaternionic charge fields Q{t',t) denotes 
the quantum spin glass order parameter field; quantum 
dynamics necessitate the solution of dynamical selfcon- 
sistency equations for the spin autocorrelation function 
< (5°°(r, t',t) > at saddle point level. Apart from the 
different interaction term in the Lagrangian, the selfcon- 
sistency equation for ri^^ resembles the LISA equation 
|]l8| for a model with nonrandom nearest neighbor hop- 
ping on a Bethe lattice with infinite coordination num- 
ber, commonly used to study strong correlations in the 
Hubbard model ||l8| . Quantum dynamics of the Q-fields 
can be important close to the T = paramagnet - SG- 
transition; yet the influence of the fermion kinetic en- 
ergy on the DoS is described sufficiently well by a static, 
replica symmetric approximation for the spin autocorre- 
lation function. For H = and half filling the Green 
function follows from 

ri = ^E^ < (ei + r,)/((e; + + H^) >$ (5) 



as given by Cy = pa,{Ei)Pa2{E2)---PaAEk), where the ^^^^ ^ j^^j^^ 



fli denote distinct replicas, do not factorize within the 
spin glass phase. For \Ek\ = E < Eg{H) we obtain 

and Beta-function B. These multi-valley correlations 
also vanish exponentially within the gap as T — > and, 
weaker than the averaged DoS though, they increase as 
either gap edge is approached at fixed low T. 
A comparison of the free energy of ordered and disordered 
phase predicts a first order transition to the paramagnetic 



The regu- 
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larized determinant 4>{H) — rrlog[l -i- ^^^^^^^2 j 
from integrating out the decoupled fermion fields, and 
H{z, y) = J{^z + VTxy). In order to proof the stabil- 
ity of the hard gap against fermionic hopping we present 
a solution of the above equation valid in the low energy 

regime. The leading contribution in a double expansion 

^2 

in Eq and frequency e; is r; = < -p- >$7sg/ giv- 
ing just back the ISGf - result. Now we note that the 
ansatz n = a ei i)satisfies eq.(||) to all orders in Eq and 
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to 0{ei) for a 



and ii) corresponds 



to a real G^ie) = ^^§t and hence to a vanishing DoS. 



Solving a set of coupled self consistency equations one 



obtains qa — I + 0[T), Xa = 



l+a 



+ 0(T), and finally 



< "p- >*„=< "p- >*/sG^= 0.568659. The last resuh 
implies that for Eqc = 2.65219 J the low energy approxi- 
mation breaks down since a diverges. On the other hand 
it is known that for Eq — 32J/(37r) the system under- 
goes a paramagnet to spin glass quantum phase transi- 
tion and has a semielliptic DoS at zero temperature , 
hence we conclude that the divergence of a at i?oc in- 
dicates the neighborhood of a metal insulator transition 
exclusively caused by random interactions. A mobility 
edge induced by either randomness in fermion hopping or 
in the interaction can interfere with the SG-based MIT. 
A transition between gapped and non-gapped insulating 
behavior appears possible too and requires further anal- 
ysis. Several different lower critical dimensions such as 
the one for Ising spin glass order at Tc > and the one 
for Anderson localization will also play a role in this con- 
text. Anderson localization of states neighboring a gap 
had been studied in Ref.( |^ ) 

Using these results for the interpretation of experiments 
one should distinguish two classes of systems. In GeCr 
, irradiated polymers [p|jl0[| , SiB ^ , ion implanted 0] 
and Czochralski grown |^ SiAs localized impurity spins 
are the main source of magnetic behavior. The interac- 
tion between them is due to direct wave function overlap, 
hence the spatial randomness in their positions translates 
itself into a broad distribution of exchange energies favor- 
ing spin glass freezing in the ground state. We believe 
that this experimental situation is well described by the 
ISGf and that the observed activated conductivity is ex- 
plained by the spin freezing induced hard gap in the DoS. 
However, in Indium-doped CdMnTe [|l| and annealed 
SiMn H] films the carrier spins interact with Manganese 
spins, may form magnetic polarons and thus create a fi- 
nite activation energy for hopping transport. In this 
mechanism is held responsible for the observed low T be- 
havior of the conductivity. Since CdMnTe is a prototype 
insulating spin glass a SG-origin of the activation energy 
can only be ruled out by checking whether the thermal 
crossover of the conductivity (in In-doped material) ap- 
pears already above the freezing temperature. If this 
were not the case an explanation of these experiments 
could be found in a spin glass of randomly interacting 
magnetic polarons. Other mechanisms discussed in the 
literature as cause for the observed activated conduction 
are a hardening of the Coulomb gap by formation of elec- 
tronic polarons 1 11 1 and nearest neighbor hopping. The 
former can be ruled out because of the observed negative 
magnetoresistance, the latter because of Efros-Sklovskii 
variable-range hopping above the crossover temperature. 



The possibility of a dimensional crossover at low T needs 
to be considered, but it would lead to a (Tq/T)^/'^ be- 
havior and not to the appearance of a hardgap. 
We acknowledge support by the SFB410 of the DFG. We 
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